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We study the box-measure correlation function of quantum states at the Anderson transition 
point with taking care of anomalously localized states (ALS). By eliminating ALS from the 
ensemble of critical wavefunctions, we confirm, for the first time, the scaling relation z(q) = 
d + 2r(g) — r(2q) for a wide range of q, where q is the order of box- measure moments and z(q) 
and r(q) are the correlation and the mass exponents, respectively. The influence of ALS to the 
calculation of z(q) is also discussed. 
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1. Introduction 

It is widely accepted that critical wavefunctions at the 
Anderson metal-insulator transition are multifractal. 1 ' 2 
Since the multifractal nature at criticality forms a ba- 
sis of the single-parameter scaling argument of the An- 
derson transition, 3 properties of critical wavefunctions 
have been extensively studied. 4 Multifractal properties 
of wavefunctions are usually characterized by the mass 
exponent r{q) or the multifractal spectrum /(a), where 
r{q) describes the system-size dependence of the local 
moments Z q of the wavefunction distribution and f(a) is 
the fractal dimension of the spatial distribution of boxes 
with the Lipschiz-Holder exponent a. It is well known 
that r(q) is related to f(a) through the Legendre trans- 
form. This relation provides a key criterion in judging 
whether numerically calculated exponents are appropri- 
ate. However, difficulty in calculating f{a) from its defi- 
nition prevents us from utilizing actually this relation for 
checking the accuracy of the exponent. For this purpose, 
we should, thus, provide other relations between expo- 
nents characterizing multifractality of a critical wave- 
function. A possible candidate is the scaling relation 
z(q) = d + 2r(q) — r(2q), where z{q) is the correlation ex- 
ponent describing the box-measure correlation function. 
The box-measure correlation function is a generalization 
of the local moment Z q and can be easily obtained in nu- 
merical calculations. Although Pracz et al. have tried to 
confirm this scaling relation for the quantum Hall tran- 
sition, 5 z(q) coincides with 2 + 2r(q) — T~(2q) only in a 
narrow range of q around q = in their result. Since the 
profile of z(q) near q = is not sensitive to details of the 
multifractal distribution, it is important to confirm the 
scaling relation for large q values. 

In order to clarify the reason why the calculated cor- 
relation exponent z(q) deviates from 2 + 2r(q) — r(2q) 
at large g's in their work, we notice the fact that some 
states are spatially "localized" due to statistical fluctua- 
tions of disorder even not in an insulating phase. These 
quantum states are called anomalously localized states 
(ALS). 6-14 Major interest in ALS is focused on them in 
a metallic phase because of the lack of theoretical tools 



for critical states. There should, however, exist ALS also 
at the critical point, because the statistical fluctuation of 
disorder at criticality easily yields a pre-localized state 
compared to ALS in the metallic regime. In fact, several 
previous works directly and indirectly support the exis- 
tence of ALS at the critical point. 15-17 Recently, we have 
shown that the ratio of the number of ALS at criticality 
to the number of the total states increases with system 
size L and saturates at a certain value for L — > oo, while 
typical states are kept to be multifractal. 18 Since the 
scaling theory and the renormalization group theory are 
based on the multifractal nature of critical wavefunctions 
(i.e., divergence of the correlation length at the transi- 
tion point), conclusions concerning properties just at the 
critical point obtained from these theories are valid only 
for typical states. Therefore, the scaling relation on the 
correlation exponent z(q) should be considered for typi- 
cal critical states. This would be a primary reason why 
the scaling relation has not been provided for large g's 
in the previous numerical work. 

In this paper, we investigate the scaling relation for 
the correlation exponent z(q) by taking into account the 
influence of ALS at the critical point. Critical wavefunc- 
tions are prepared in two-dimensional symplectic sys- 
tems which are described by the SU(2) model. 19 In order 
to obtain z(q) for typical critical states, we compose an 
ensemble of critical wavefunctions from which ALS are 
eliminated. The exponent z(q) calculated for such an en- 
semble satisfies the scaling relation for large q values. 
We quantify how well the ensemble from which ALS are 
eliminated describes statistical properties of typical crit- 
ical wavefunctions by a parameter V* . It is found that 
the range of q for which z(q) satisfies the scaling rela- 
tion depends on T* . This paper is organized as follows. 
In §2, we give a quantitative definition of ALS at the 
critical point based on the idea that ALS at criticality 
do not show multifractality. In §3, the scaling relation 
for the correlation exponent is reminded in the scaling 
argument for the box-measure correlation function. We 
briefly explain, in §4, the SU(2) model and a numerical 
method to obtain eigenstates of this model. Results of 
our numerical study are given in §5. The insensitive na- 
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ture of z(q) near q — to ALS is also discussed in this 
section. Section 6 is devoted to conclusions. 

2. Definition of Anomalously Localized States 

In order to study the influence of ALS quantitatively, 
we employ a definition of ALS proposed in ref. 18. This 
definition is based on the idea that ALS are not mul- 
tifractal as a consequence of their localized nature. At 
first, we introduce the box- measure correlation function 
G q (l,L,r) defined by 20 - 21 

r b b r 

where p m = J2ieb(i) IV'il 2 and Vb r (l) = J2teb r (i) ^ are 
the box measures for wavefunction amplitudes ipi, in a 
box b(l) of size I and in a box b r (l) of size I fixed distance 
r — I away from the box b(l), respectively. Nb (or Nf, r ) is 
the number of boxes 6(7) [or b r (l)], and the summation 
J2b ( or J2b ) i s taken over all boxes b(l) [or b r (l)} in 
the system of size L. If a wavefunction is multifractal, 
G q (l, L, r) should behave as 21 

G g (l,L,r) ccl x{q ^L- v{q \- z{q \ (2) 

where x(q), y(q), and z(q) are exponents describing mul- 
tifractal correlations of the amplitude distribution. This 
relation is sensitive to ALS as demonstrated in ref. 18 
and then suitable for defining ALS. To find the / and r 
dependences of G q (l, i, r), we concentrate on the follow- 
ing functions, 

Q q (l) =G q (l,L,r = l) «/*<')-*('), (3) 

and 

R q (r)=G g (l = l,L,r)ocr- z ^. (4) 

In order to quantify non-multifractality of a specific 
wavefunction, it is convenient to introduce variances 
Var(logQ 2 ) and Var(logi? 2 ) from the linear functions of 
log? and logr, logQ 2 (0 = [x(2) - z(2)]log/ + c Q and 
log i? 2 (r) = —2(2) logr + cr, respectively, calculated by 
the least-square fit. From these variances, a quantity F 
is defined by 

r(L, A) - AVar(logQ 2 ) + Var(lo gj R 2 ), (5) 

where A is a factor to compensate the difference between 
average values of Var(logQ 2 ) and Var(logi? 2 ). Using T 
given by eq. (5), the quantitative and expediential defi- 
nition of ALS at criticality is presented by 

F > T*, (6) 

where F* is a criterial value of T to distinguish ALS from 
multifractal states. We can compose a refined ensemble 
by eliminating ALS from a set of critical wavefunctions. 
The quality of the refined ensemble is controlled by T* . 

3. Scaling Relation for the Correlation Expo- 
nent 

In this section, we give a brief explanation of the scal- 
ing relations between the exponents x(q), y(q), z(q), and 
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r(q). The mass exponent r(q) which is commonly used 
in the multifractal analysis is defined by 

^(o^E^o Kr(9) - ( ? ) 

b 

The local moment Z q (l) is related to the correlation 
function G q (l 7 L,r). Comparing definitions of G q and Z q 
[eqs. (I) and (7)], we have 

Gq (l,L,r = l) = ^±. (8) 

Therefore, the relation N b = {L/l) d and eqs. (2) and (7) 
lead 

x(q)-z(q)=y(q)=d + T{2q). (9) 
For the case of r = L, eq. (1) gives 



G q (l,L,r = 




Since G q (l, L,r = L)(x from eq. (2), the 

exponents are related as 

x(q)=y(q) + z(q) = 2d + 2r(q). (11) 

From eqs. (9) and (11), we obtain 

x(q) = 2d + 2 T (q), (12) 

y(q) = d + r(2q), (13) 

and for the correlation exponent 

z{q) = d + 2r{q)-T{2q). (14) 

The last relation with q = 1 is equivalent to the well- 
know scaling relation 77 = d — Z? 2 , 22 where 77 is the ex- 
ponent describing the two-particle correlation function 
S(r) = p(\ip(0)\ 2 \ip(r)\ 2 ) , p is the density of states, and 
D2 is the correlation dimension of critical wavefunctions. 
This equivalency can be understood from the relations 
z(l) = r/ because of G\(l = l,L,r) oc S(r), r(l) = 0, 
and t(2) = D 2 . 

The above relations (12)-(14) hold for any q in prin- 
ciple if the wavefunction is multifractal. However, the 
previous work confirming the scaling relation eq. (14) 
for the quantum Hall transition 5 shows that the correla- 
tion exponent z(q) coincides with 2 + 2r(q) — r(2q) for 
a narrow range of q (\q\ < 1.5). As mentioned in §1, we 
consider that the disagreement at large |g|'s is a conse- 
quence of ALS. It should be noted that the agreement 
for small |g|'s is rather trivial. Even if the wavefunc- 
tion is not multifractal, the exponents satisfy the scal- 
ing relation eqs. (12)-(14) for q = 0, because r(0) = — d 
and x(0) = j/(0) = z(0) = for any Pb(i)- Increasing 
\q\ slightly from zero, the correlation exponent z(q), for 
example, is well approximated by the parabolic form 

z(q) = 2(a Q -d)q 2 , (15) 

where «o is the Lipschiz- Holder exponent giving the max- 
imum of f(a). If a is insensitive to details of the distri- 
bution of measures, z(q) approximately satisfies eq. (14) 
near q = 0, because r(q) is given by (1 — q) [q(a — d) — d] 
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in the parabolic approximation. We should, thus, con- 
firm the scaling relations for large \q\'s. Since the expo- 
nent z(q) characterizes directly the spatial correlation of 
wavefunction amplitudes, we focus only on the relation 
(14) hereafter. 

4. System and Numerical Method 

Considering the advantage of system sizes, we fo- 
cus our attention on the Anderson transition in two- 
dimensional electron systems with strong spin-orbit in- 
teractions, in which systems have no spin-rotational sym- 
metry but have the time-reversal one. Hamiltonians de- 
scribing these systems belong to the symplectic ensem- 
ble. Among several models belonging to this universality 
class, we adopt the SU(2) model because of its small scal- 
ing corrections. The Hamiltonian of the SU(2) model 19 
is given by 



H 



E 



C Ci 



(16) 



t,3 



where cj (cj) is the creation (annihilation) operator act- 
ing on a quaternion state vector, R^ is the quaternion- 
real hopping matrix clement between the sites i and j, 
and £$ denotes the on-site random potential distributed 
uniformly in the interval [—W/2,W/2]. (Bold symbols 
represent quaternion-real quantities.) The matrix ele- 
ment Rij is given by 

Rij — cos a.i j cos Pi j t° + sin 7^- sin r 1 

— cos 7y sin f3ij t 2 + sin a,j cos (3ijT 3 , (1 7) 

for the nearest neighbor sites i and j, and R^ = for 
otherwise. Here, t m (ji = 0,1,2,3) is the primitive ele- 
ment of quaternions. 23 Random quantities cty and 7^ 
are distributed uniformly in the range of [0, 2n), and 
(3ij is distributed according to the probability density 
P(/3)dP = sin(2/3)d/3 for < (3 < tt/2. Randomly distri- 
buted hopping matrix elements shorten the spin relax- 
ation length which is a dominant irrelevant length scale, 
and scaling corrections can be negligible in the SU(2) 
model. It is known that the localization length exponent 
v of this model is 2.73 ±0.02 and the critical disorder W c 
is 5.952V at E= l.OV. 19 

Critical wavefunctions of the SU(2) model have 
been calculated by using the forced oscillator method 
(FOM) 24 extended to the eigenvalue problem of 
quaternion-real matrices. Of course, the Hamiltonian 
eq. (16) can be represented by complex numbers, and 
we can use the usual FOM for complex Hcrmitian matri- 
ces to solve the eigenvalue problem. The modified FOM 
for quaternion-real matrices, however, enables us to cal- 
culate eigenvalues and eigenvectors within about a half 
of CPU time. 25 It should be remarked that the obtained 
eigenvector is a quaternion-real vector. This vector rep- 
resents two physical states simultaneously, which corre- 
spond to the Kramers doublet. Since the amplitude dis- 
tribution of these degenerate states are the same, we an- 
alyze one of the calculated Kramers doublet. 

5. Results 

In order to study the correlation exponent z(q) for crit- 
ical wavefunctions of the SU(2) model, we calculate 10 4 
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Fig. 1. Calculated correlation exponent z(q) (open circles) and 
the right-hand side of cq. (14) [z T (q) = d + 2r(q) - r(2q)] (filled 
circles). Figures labelled by (a) and (b) correspond to the results 
for the original ensemble and the refined ensemble, respectively. 
The dashed lines show the parabolic approximation of z(q) given 
by eq. (15) with qq = 2.18. 



critical wavefunctions of this model at E = l.OV and 
W = 5.952V by the FOM. Each eigenstate is obtained 
for a single disorder realization. Periodic boundary con- 
ditions are imposed in the x and y directions in systems 
of size L = 120. 

Figure 1(a) shows the q dependence of the correlation 
exponent z{q) and the right-hand side of eq. (14), i.e., 
d+2r{q) — r(2q)[= z T (q)], calculated for the whole set of 
the prepared critical wavefunctions (the original ensem- 
ble). The exponents z(q) and r(g) are calculated from the 
geometric means of R q (r) and Z q (l) by the least-square 
fit, respectively. It is found that z(q) deviates from z T (q) 
at large l^l's. The degree of agreement is similar to the 
result of the previous work for the quantum Hall transi- 
tion. 5 We see that z(q) for q > 5 exceeds the value of 2 
which is the theoretical bound of z(q). 

It should be noted that the scaling relation eq. (14) 
is valid for typical critical wavefunctions. As shown in 
our previous work, 18 the distribution function of T de- 
fined by eq. (5) has its peak at T = 0, which implies 
that typical critical wavefunctions are multifractal. Since 
the geometric mean usually represents the typical value 
(mode value), it seems that the geometric mean (R q (r)) 
used in our calculations gives the correct exponent z(q) 
of multifractal critical wavefunctions without influence 
from ALS. However, the distribution function of T de- 



4 J. Phys. Soc. Jpn. 



Full Paper 



Author Name 



fined for the correlation function G q becomes broad with 
increasing \q\, and then the geometric mean could not 
represent the typical value for an ensemble composed of 
an insufficient number of elements due to enhanced sta- 
tistical fluctuations. In fact, we have confirmed that the 
correlation exponent calculated from the typical profile 
of R q (r) largely deviates from the exponent computed 
from (R q (r)) for large \q\. Although the correlation ex- 
ponent appeared in eq. (14) should be computed from the 
typical profile of R q {r) [R q yp (r)] : it is actually quite dif- 
ficult in numerical calculations to calculate precise value 
of z(q) from Rf p (r). 

In this paper, we show that the elimination of ALS 
from the original ensemble of critical wavefunctions is 
efficient technique to obtain accurate value of z(q) even 
for large \q\. At first, we compose a refined ensemble by 
eliminating ALS from the original ensemble. The refined 
ensemble is characterized by the parameters A = 3 and 
T* = 0.06, and contains 6962 critical wavefunctions. The 
correlation exponent calculated from (R q (r)) for this re- 
fined ensemble is presented in Fig. 1(b). The exponent 
z(q) is in agreement with z T (q) for a wider range of q, and 
does not exceed the theoretical bound 2. This demon- 
strates the efficiency of the elimination scheme of ALS to 
calculate precise exponents defined at the critical point. 
We find, however, that the exponent z(q) deviates from 
z T (q) for q < — 2 even for the refined ensemble. What 
is even worse, z{q) becomes negative for q <C —2. To 
see the reason of this unreasonable result, let us con- 
sider a simple case of R q (r), namely, R q (r) for quantities 
Hi distributed uniformly in the range of [0, 1] without 
any spatial correlations. From the definition of R q (r) [or 
G g (i,L,r)], we have 

r ^ = wnX £ ^ (18) 



where N is the total site number, Y^jerU) represents the 
summation over sites j away from the site i by a distance 
r, and N r is the number of such sites. The uncorrelated 
distribution of makes it possible to replace /x? by its 
average value [i q , and we obtain 



R q (r) 



m 2 

NN r 



E E = ^f- 

i jer(i) 



(19) 



From this relation, it seems that the exponent z(q) be- 
comes zero for any q, because ~jjfl does not depend on 7-. 
This is, however, not true for negative q. The distribution 
function of [i q is given by Q{t)l\qt^-^li\ for positive q 
and [1— 9(t)]/[\q\t( q ~ 1 ^ q ] for negative q, where t = [i q and 
8(t) is the step function defined as 6(t) = 1 for < t < 1 
and 6(t) = for t > 1. For q > 0, the distribution is 
truncated at t = 1 and exponent (q — l)/q in the dis- 
tribution function is less than unity, which gives a finite 
average value of \x q . On the contrary, the distribution 
function for negative q extends over the range of [1, oo). 
The exponent (q — l)/q becomes less than 2 for q < —1, 
which leads an infinite average value of fi q . Thus we can- 
not define the exponent z(q) from eq. (19) for q < — 1. 
The divergence of ~jjfl is a consequence that the distribu- 
tion function of /j is finite at /i = 0. Since a distribution 




Fig. 2. Multifractal spectra /(«) for the original ensemble (cir- 
cles) and the refined ensemble (crosses). Crossed error bars in- 
dicate representative values of errors in /(«) and a. The dashed 
line indicates the parabolic approximation eq. (20) with ao = 
2.18. 



of multifractal measures is not uniform (but log-normal 
approximately) and spatial correlations of them must be 
taken into account, the above argument cannot be ap- 
plied straightforwardly to refined ensembles of critical 
wavefunctions. Nevertheless, the log-normal distribution 
function taking relatively large values near [i = induces 
a numerical instability for q < 0. This instability leads 
the unreasonable result for q <C — 2. 

We see from Fig. 1 that the exponent z T (q) is not sen- 
sitive to ALS. This is consistent with the fact that the 
local moment Z q (l) is not enough to distinguish ALS 
from multifractal wavefunctions. The multifractal spec- 
tra calculated for the original and the refined ensembles 
also indicate this fact as shown in Fig. 2, where /(a) has 
been calculated by the g-microscope. 26 Especially, these 
two spectra almost collapse in the vicinity of ao giving 
the maximum of f(a), which corresponds to q — 0. Since 
the profile of f{a) close to a = ao can be approximated 
by the parabolic form 



/(<*) = 2 - 



(op - a) 2 
4(a - 2) 



(20) 



this implies that ao is hardly affected by ALS. (Dashed 
line in Fig. 2 shows this parabolic approximation with 
a = 2.18.) Using this value of ao, we can draw the 
parabolic approximation of z(q) given by eq. (15) as 
shown in Fig. 1. The insensitivity of ao to ALS leads 
insensitive z(q) to ALS near q = 0. In fact, z(q) shown 
in Fig. 1 coincides with z T (q) and the parabolic approx- 
imation in the vicinity of q = even for the original 
ensemble. 

Apparently, the quality of the refined ensemble de- 
pends on the choice of the value of T* in eq. (6). Figure 3 
exhibits this quantitatively. The longitudinal axis q* of 
this figure represents the value of q (> 0) above which 
\z(q) — z T (q)\ is larger than 0.04. The data points indi- 
cated by arrows correspond to the original ensemble and 
the refined ensemble used in Figs. 1 and 2. This result 
indicates that the exponent z(q) coincides with z T (q) in 
the wider range of q when the degree of refining the en- 
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Fig. 3. Quantity q* representing the quality of the refined ensem- 
ble as a function of T*. q* is defined as \z(q) — z T (q)\ > 0.04 for 
q > q* . Arrows indicate the data points corresponding to the 
refined ensemble used in Fig. 1 ( or 2) and the original one. 

semble becomes higher. Namely, the high quality refined 
ensemble ensures the scaling relation eq. (16) with large 
q. It seems that q* becomes infinite when T* goes to zero. 
The quantity q* , however, does not diverge actually. Too 
small r* makes the quality of the refined ensemble rather 
worse, because the number of samples in the refined en- 
semble becomes insufficient to provide a good statistics. 

6. Conclusions 

We have investigated the box-measure correlation 
function G q (l, L, r) of critical wavefunctions of the SU(2) 
model which belongs to the two-dimensional symplectic 
class and confirmed the scaling relation for the corre- 
lation exponent z(q) in a wide range of q. It is found 
that ALS prevent us from calculating accurate values of 
the correlation exponent with large \q\'s. To obtain z(q), 
we eliminated ALS from the original ensemble of criti- 
cal wavefunctions and composed a refined ensemble for 
which the geometric mean of the box-measure correla- 
tion functions represents properly the typical profile of 
G q (l,L,r). It is also found that the range of q within 
which the calculated z(q) satisfies the scaling relation 
becomes wider when the quality of the refined ensemble 
becomes better (namely, decreasing T*). We have elu- 
cidated that a direct calculation of the exponent z(q) 
encounters a numerical instability for q <C —2, which 
originates from the fact that the distribution function of 
multifractal measure [i has relatively large values near 
/i = 0. Some ingenious techniques in numerical calcula- 
tions would be required to obtain precise values of z(q) at 
negative q's. In the present paper, the dependence of the 
correlation function G q (l,L,r) on the box-box distance 
r has been examined only for I = 1. If we can deal with 
huge systems, G q (l, L, r) with I > 1 would provide a more 
precise value of z(q) for q <C —2, because the numerical 
instability at negative q's might be moderated due to an 
average effect within each box. In our demonstrations by 
using the SU(2) model, the value of q to eliminate ALS 
has been fixed at 2, while G q has been analyzed for many 
q values. Essentially, the choice of the value of the mo- 



ment q for the elimination does not influence the quality 
of the obtained refined ensembles. It is, however, better 
to use a positive small value of q for ensuring a numerical 
stability. 

The composition of the refined ensemble is efficient 
also for evaluating other properties at the critical point 
than the scaling relation treated in this paper. We em- 
phasize that the scaling relation eq. (14) can be used for 
judging how well the refined ensemble represents typi- 
cal critical wavefunctions. It is also interesting to study 
statistical properties of the ensemble of ALS, namely, 
the complementary set of the refined ensemble. Further 
investigations of ALS provide us deep insight into the 
nature of criticality. 
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